Quantum criticality of a Bose gas in an optical lattice near the Mott transition 
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We derive the equation of state of bosons in an optical lattice in the framework of the Bose- 
Hubbard model. Near the density-driven Mott transition, the expression of the pressure P{n, T) 
versus chemical potential and temperature is similar to that of a dilute Bose gas but with renor- 
malized mass m* and scattering length a* . m* is the mass of the elementary excitations at the 
quantum critical point governing the transition from the superfiuid phase to the Mott insulating 
phase, while a* is related to their effective interaction at low energy. We use a nonperturbative 
renormalization-group approach to compute these parameters as a function of the ratio t/U be- 
tween hopping amplitude and on-site repulsion. 
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Introduction. Quantum phase transitions play a cru- 
cial role in many strongly-correlated systems ranging 
from quantum antiferromagnets to heavy-fermion mate- 
rials and high-T c superconductors [ ] . Although these ' •■? 
transitions occur at zero temperature, theory predicts 
that the finite-temperature thermodynamics in the vicin- 
ity of a quantum critical point (QCP) is described by 
universal scaling relationships up to rather high temper- 
atures [2, 3]. While understanding quantum criticality 
in strongly-correlated systems is often a challenge, both 
experimentally and theoretically, cold atoms offer clean 
systems for a quantitative and precise study of quan- 
tum phase transitions. Quantum criticality in cold atoms 
has attracted increasing theoretical interest in the last 
years [ ] , and the experimental observation of the quan- 
tum critical behavior of a two-dimensional Bose gas in 
an optical lattice near the vacuum-supcrfluid transition 
has recently been reported [5]. 

In ultracold gases, strong correlations can be achieved 
by tuning the atom-atom interactions by means of a Fes- 
hbach resonance, or by loading the atoms into an optical 
lattice [ ]• In the latter case, by varying the strength of 
the lattice potential, it is possible to induce a quantum 
phase transition between superfiuid and Mott insulating 
ground states in a Bose gas [ ]. The main features of 
this transition can be understood in the framework of 
the Bose-Hubbard model, which describes bosons hop- 
ping on a lattice with an on-site repulsive interaction [ ] . 

The density-driven Mott transition belongs to the same 
universality class as the transition between the vacuum 
and the superfiuid phase in a dilute Bose gas: both tran- 
sitions are governed by the same QCP [3, 8]. It follows 
that the pressure must take the same (universal) form 
near the QCP (up to some nonuniversal parameters as 
will be explained in detail below). In this Letter, we 
prove this result by computing the pressure P{n, T) (as 
a function of chemical potential and temperature) in vari- 
ous limits using a nonperturbative renormalization-group 
approach to the Bose-Hubbard model [9, 10]. 
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FIG. 1. (Color online) Phase diagram of the Bose-Hubbard 
model on a cubic lattice showing the Mott insulators (MI) 
with density n — (vacuum) and n = 1, as well as the sur- 
rounding superfiuid phase (SF). Point C at the tip of the Mott 
lobe shows the multicritical point where the transition occurs 
at fixed density n = 1. Away from this point, the transition 
is driven by a density change. The zero-temperature pressure 
P{n, T = 0) along the dotted line AB is shown in Fig. 4. The 
finite-temperature pressure P(// C ,T) at point B is shown in 
Fig. 5. (In all figures of the Letter, we take U as the energy 
unit.) 



The equation of state of a dilute Bose gas is usually de- 
rived from a low-density expansion. To next-to-leading 
order in the parameter mo?\i, the pressure P(p, T) de- 
pends only on the boson mass m and the s-wave scatter- 
ing length a (and not on the microscopic details of atom- 
atom interactions). The equation of state can also be un- 
derstood from the point of view of phase transitions [3] . 
By varying the chemical potential from negative to pos- 
itive values at zero temperature, one induces a quantum 
phase transition between a state with no particles (vac- 
uum) and a superfiuid state. Elementary excitations at 
the QCP are free bosons of mass m and the interaction 
between two elementary excitations is determined by the 
scattering length a in the low-energy limit. Thus the 
equation of state P(/i, T) is uniquely determined by the 
QCP which therefore controls the thermodynamics of the 
dilute Bose gas. 



The advantage of this point of view is that it allows us 
to understand the density-driven Mott transition in the 
Bose-Hubbard model along similar lines, using the fact 
that it belongs to the dilute Bose gas universality class. 
At the QCP between the Mott insulator and the super- 
fluid phase, the elementary excitations are quasi-particlcs 
with effective mass m* and their mutual interaction is de- 
scribed by an effective scattering length a*. Thus when 
the boson density n slightly differs from the commensu- 
rate density n c of the Mott insulator, the excess density 
\n — n c \ of particles (or holes) behaves as a dilute gas of 
quasi-particles [8] (with effective parameters m* and a*) 
with an equation of state P(/z, T) which must be the same 
as the one of the dilute Bose gas up to some nonuniversal 
parameters such as m* and a*. 

Our main results can be summarized as follows. Near 
the density-driven Mott transition and for temperature T 
smaller than the hopping amplitude t, the pressure takes 
the form 

P((i, T) = P C + n c 5[i + m* 3/2 T 5 ^P (±^, ±m*a* 2 S^ , 

(1) 
where P c and n c are the pressure and density (mean num- 
ber of bosons per site) at the QCP \x = \i c , respectively, 
and Sfj, = j-i — fi c . P(x, y) is a universal scaling function 
characteristic of the (three-dimensional) dilute Bose gas 
universality class. The + (— ) sign in (f) corresponds to 
particle (hole) doping. Equation (f ) is valid everywhere 
near the superfluid-Mott-insulator transition except in 
the close vicinity of the multicritical points where the 
transition takes place at fixed (commensurate) density 
(Fig. f). For fi c = P c = n c = 0, m* = m and a* = a, 
Eq. (1) reproduces the low-density expansion of a dilute 
gas of bosons with mass m and scattering length a. In 
this Letter, we discuss the equation of state (f ) in two 
different limits. 

In the zero-temperature superfluid phase, 
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P(jji,T = 0) = P c + n c Sn 
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where the ellipses denote higher-order terms in the ex- 
pansion parameter m*a* 2 |<5«|. Taking the derivative of 
Eq. (2) wrt /x, we obtain 

m*Sfi ( 16 
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where n = dP/dfj, is the boson density. Equations (2) 
and (3) are similar to the known results for a dilute 
Bose gas. The correction to the "mean-field" result 
Pq = m*(Sfi) 2 /8ira* is known as the Lee-Huang- Yang 
correction [ ]. 

At finite temperature and for /i = « c (quantum critical 
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(4) 



for T < t, where t is the hopping amplitude between 
neighboring sites [See Eq. (5) below]. Once T is of order 
of the boson dispersion t, the lattice starts to play a role 
and the universal description in terms of quasi-particles 
(with parameters m* and a*) breaks down. Equation (4) 
agrees with the well-known expression of the pressure in 
the dilute Bose gas when |/z| <C T [3]. 

Equation of state. To derive Eqs. (2) and (4), we start 
from the action of the Bose-Hubbard model, 



S = l drlY\r T {d T - ti)i>r 
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where ^) r {j) is a complex field and r € [0, /3] an imaginary 
time with f3 = 1/T the inverse temperature, {r} denotes 
the N sites of the lattice, U the on-site repulsion, and 
t the hopping amplitude between nearest-neighbor sites 
(r,r'). We set h = ks = 1 and take the lattice spacing 
as the unit length throughout the Letter. 

The NPRG approach [9, 10] allows us to compute the 
effective action T [</>*,</>] defined as the Legendre trans- 
form of the thermodynamic potential — In Z[J* , J]. J r is 
an external source which couples linearly to the bosonic 
fielder and0 r (r) = 8lnZ[J*, J]/5J?(t) = (M T )) is the 
superfluid order parameter. Thermodynamic properties 
of the system can be derived from the effective poten- 
tial defined by V(n) = (0N)- 1 T[())*,<p] with n = \(j>\ 2 
and <j) a constant (uniform and time-independent) field. 
Its minimum determines the condensate density no an d 
the pressure P = — V(n ) in the equilibrium state. The 
single-particle propagator G = — r' 2 ^ 1 is related to the 
two-point vertex r' 2 ' defined as the second-order func- 
tional derivative of T. 

At the QCP between the superfluid phase and the Mott 
insulator, the effective action takes the form 



rr«f,, 



dr / d 6 r 
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(6) 



(with j3 — > oo) up to a constant (field-independent) term. 
Since we are interested in the low-energy limit, we con- 
sider the continuum limit where r becomes a continuous 
variable. The ellipses denote higher-order (in derivative 
or field) terms. Equation (6) is valid at a generic QCP 
where the transition is driven by a density change. At a 
multicritical point, where the transition occurs at fixed 
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FIG. 2. (Color online) Zc, A and Z 4 vs t/U at the QCP 
between the superfluid phase and the Mott insulator n — 1. 
The + and — signs refer to the upper and lower parts of the 
transition line. 
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FIG. 3. (Color online) Effective mass m* and scattering 
length a* vs t/U at the QCP between the superfluid phase 
and the Mott insulator n = 1. The solid (red) line shows the 
scattering length a — [8iv(t/U + A)] -1 (see text). 



(commensurate) density, Zc vanishes and one should ex- 
plicitly include a d 2 term; the quantum phase transition 
is then in the university class of the 0(2) model [9, 10]. 
From Eq. (6), we can identify the elementary excita- 
tions at the QCP. Zc is negative on the lower part of the 
transition line (for a given Mott lobe). In that case, it 
is convenient to perform a particle-hole transformation 
4> <H- <fi* (which changes the sign of the d T term in (6)). 
We can then define a quasi-particle field cf> = \Zc\ 1 ^ 2 4> 
and rewrite the effective action as 



r$*,$\ = 



(It 
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Za 



4^' 



(7) 



(8) 



with m = l/2t the effective mass of the free bosons mov- 
ing on the lattice. The elementary excitations are quasi- 
particles with a quadratic dispersion law and mass m*. 
They are particle-like if Zc > and hole-like if Zc < 0, 
and the quasi-particle weight is \Z C \~ X ■ The interaction 
between two quasi-particles is determined by the effective 
scattering length a*. 



0.16 



0.12 - 



0.08 



0.04 




-0.1 -0.05 



0.05 0.1 



0.15 




0.0015 



0.001 



FIG. 4. (Color online) Pressure P vs chemical potential \i 
along the dotted line in Fig. I, as obtained from the NPRG 
approach (the inset shows the density n — dP/dfi). The bot- 
tom figures show the behavior near the Mott insulating phases 
n — and n — 1. The dashed (green) line corresponds to 
Eq. (2) and the dash-dotted (blue) one to the "mean-field" 
result P = P c + n c S/i + (Sfi 2 )m* /8ira* . 



At the QCP between the superfluid phase and the 
vacuum (the Mott insulator with vanishing density), 
Zc = Za = 1 (the single-particle propagator is not 
renormalized [l 11 ]) so that m* = m = l/2t. Further- 
more, the interaction constant A = 87rta can be calcu- 
lated analytically and related to the scattering length 
a = [8ir(t/U + A)}- 1 (A ~ 0.1264) of the bosons moving 
on the lattice [10], which gives a* — a. For a generic 
QCP between the superfluid phase and a Mott phase 
with nonzero density (h c = 1, 2, • • • ), the determination 
of to* and a* is much more difficult as it requires to solve 
a many-body problem. Zc, Za and A can be obtained 
from the numerical solution of the NPRG equations [ ]. 
Figures 2 and 3 show Zc, Za, A, 777* and o* as a function 
of t/U for the QCP separating the superfluid phase from 
the Mott insulator n = 1. The vanishing of Zc at the 
multicritical point implies that 777* vanishes and a* di- 
verges when we approach the tip of the Mott lobe (point 
C in Fig. 1). Note that it is numerically difficult to deter- 
mine 777* and a* for t/U < 0.015 due to the degeneracy 
between the states with density n = p and 77 = p + 1 
when fi/U — p (p integer) and t = 0. 

We are now in a position to verify that near the QCP 
the zero-temperature pressure in the superfluid phase is 
given by Eq. (2) [12]. Note that higher-order terms ne- 
glected in Eq. (6) do not contribute to the pressure to 
the order considered [13]. Figure 4 shows P(/j,,T — 0) at 
fixed t/U and for a density varying between and 1 (see 
the dotted line AB in Fig. 1). Near the Mott insulating 
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FIG. 5. (Color online) Pressure P(/i c ,T) vs temperature 
T [t — 0.02(7]. The dash-dotted (blue) line corresponds to 
Eq. (4). The inset shows a log-log plot and the T 5 ' 2 depen- 
dence of P(/i c ,T) - P c . 



P(n, T — 0) [Eq. (2)] could be observed. Recently, quan- 
tum criticality of a two-dimensional Bose gas near the 
vacuum-superfiuid quantum phase transition has been 
observed at finite temperature, thus explicitly demon- 
strating that the equation of state of cold atomic gases 
gives direct information about the QCP [5]. We expect 
that similar measurements will allow to observe quan- 
tum criticality near the (nontrivial) QCP governing the 
superfluid-Mott-insulator transition, and verify our pre- 
dictions for the effective mass m* and scattering length 
a* in three-dimensional Bose gases. 

We would like to thank X. Leyronas, F. Chevy and C. 
Salomon for discussions. 



phases n — and n = 1, we find a very good agreement 
between the NPRG result and Eq. (2). As we approach 
the tip of the Mott lobe located at (£ C ,/2 C ), m*a* 2 di- 
verges and the domain of validity \Sfj,\ <C l/m*a* 2 of 
Eq. (2) shrinks to zero. For t = i c and fi ~ /2 C , the pres- 
sure varies as P(n, T = 0) ~ P c + n c 5[i + a([i — fl c ) 4 with 
a a constant. 

The finite-temperature pressure P(fJ. C i T) at point B in 
Fig. 1 is shown in Fig. 5. We obtain a very good agree- 
ment with Eq. (4) as long as lattice effects can be ignored, 
i.e. T < t. One can also compute the temperature de- 
pendence of the pressure for /i ^ fi c . In all cases, one 
recovers the known results of the dilute Bose gas [3] but 
with m and a replaced by m* and a*. 

Conclusion. We have shown that the pressure P(/k, T) 
of a Bose gas near the density-driven Mott transition 
takes a universal form, analog to that of a dilute Bose gas, 
once a certain set of low-energy parameters have been 
fixed. These parameters characterize the elementary ex- 
citations at the QCP and their interactions, and can be 
systematically computed using the NPRG approach. In 
three-dimensions (i.e. above the upper critical dimen- 
sion d c = 2 of the T = superfluid-Mott-insulator tran- 
sition), this set is formally infinite but to leading order 
and for T <C t only the mass m* and scattering length a* 
of the elementary excitations is needed [Eq. (1)]. At the 
upper critical dimension d c — 2, our analysis still holds; 
universality is even stronger, as only m* and a* need to 
be known (to all orders in the small parameter m*a* 2 \8fi\ 
and for T <C i) to determine the equation of state. 

Recent experiments in cold atom gases have shown that 
it is now possible to deduce the pressure P(n, T) of an 
homogeneous infinite system from the doubly-integrated 
in situ density profile J dxdyn(x,y,z) of an harmoni- 
cally trapped gas [ ]. The zero-temperature equation 
of state (including the Lee-Huang- Yang correction) of 
a homogeneous Bose gas of 7 Li atoms has been mea- 
sured using this technique [15]. Experiments with atoms 
loaded in an optical lattice [ ] are now approaching 
the low-temperature regime where our predictions for 
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